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Abstract 

We invoke universal Chern-Simons theory to analytically calculate the exact free en¬ 
ergy of the rehned topological string on the resolved conifold. In the unrehned limit 
we reproduce non-perturbative corrections for the resolved conifold found elsewhere in 
the literature, thereby providing strong evidence that the Chern-Simons / topological 
string duality is exact, and in particular holds at arbitrary N as well. In the rehned 
case, the non-perturbative corrections we hnd are novel and appear to be non-triviah 
We show that non-perturbatively special treatment is needed for rational valued de¬ 
formation parameter. Above results are also extend to rehned Chern-Simons with 
orthogonal groups. 
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1 Introduction 


By now, the large N duality between SU{N) Chern-Simons on and the (pertur¬ 
bative) topological string on the resolved conifold of [1] is well established. Over the 
years this duality has been extended in several ways, like for example to SO/Sp gauge 
groups [2], leading at large N to orientifolds of the topological string. 

More recently, triggered by work of Nekrasov on A/" = 2 supersymmetric gauge 
theory [3], it became clear that there should exist a sort of rehned topological string, 
being a one-parameter deformation of the usual topological string [U |5] (for a brief 
exposition, see 0)- The deformation parameter is usually denoted as f3. In turn, this 
led Aganagic-Shakirov miH] to propose a rehned Chern-Simons theory yielding at large 
N the free energy of the rehned topological string on the resolved conifold. Of course, 
at = 1 the original large N dualities are recovered. In fact, the rehned Chern-Simons 
theory can be dehned for all ADE groups [9]. In particular, for this leads at large 
to a rehnement of topological string orientifolds. 

Perhaps less known, in a series of works a novel universal formulation of Chern- 
Simons theory on has been put forward dQi [m [12]. Here, universal refers to 
the fact that all the partition functions of Chern-Simons with arbitrary classical or 
exceptional simple gauge groups can be recovered from the universal Chern-Simons 
partition function under specialization of parameters. Quite surprisingly, the universal 
Chern-Simons theory does not only include the usual Chern-Simons theories, but, after 
some extension of range of parameters, as well the rehned versions thereof, as shown 
in [13] . 

By construction, the universal Chern-Simons partition function constitutes an in¬ 
tegral representation of the partition functions and thereby provides an analytic con¬ 
tinuation in the parameters, e.g., simple Lie groups are now parametrized by the two- 
dimensional Vogel’s plane [H]. This generalizes the old N —)■ —N duality of gauge 
theories with orthogonal and symplectic groups [15], and leads in particular to a sug¬ 
gestion of an extension of gauge/string duality to exceptional groups [IB]. Furthermore, 
the integral representation is very well suited to study non-perturbative aspects of the 
large N duality, see [niEB]. 

In this work we continue this line of research by using the universal Chern-Simons 
integral representation to analytically calculate non-perturbative corrections to the 
rehned Chern-Simons theory at large N, thereby proposing the non-perturbative com- 
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pletion of the refined topological string on the resolved conifold. However, one shonld 
keep in mind that for integer N the non-perturbative corrections vanish, as the analyti¬ 
cally continued partition function at integer N coincides with the initial Chern-Simons 
partition function. At /9 = 1 we can compare to recently renewed efforts to hnd the 
non-perturbative completion of the topological string, in particular to [ni[l8]. How¬ 
ever, we like to stress that our approach is entirely analytic and does not, in contrast 
to other works in the literature, rely on any subtle combination of quantization, ap¬ 
proximation and numerics. In fact, we do not even perform a genus (or more generally 
trans-series) nor large N expansion, but directly recover for SU{N) gauge group the 
refined Gopakumar-Vafa expansion of the resolved conifold (similar as previously in 
for the unrehned case) and, most importantly, a non-perturbative completion 
thereof, via simple residue calculations. In particular, one may see this as the simplest 
example of a true {i.e., independent of N) gauge / string duality. Similar results are 
obtained for SO{N) gauge group, leading to the prediction of the non-perturbative 
completion of orientifolds of the rehned topological string on the resolved conifold. 

The outline is as follows. In the next section we will recall the basic definitions 
of universal Chern-Simons theory. In particular, we will rewrite the known integral 
representation in a more convenient form, see section [2711 which will allow us to express 
the integral representation directly as a sum of residue, or, alternatively in terms 
of multiple sine functions, as discussed in section 12.21 In the following two sections 
examples are discussed. Namely, in section E] we calculate the universal partition 
function for unitary groups as a sum of residue, yielding the rehned Gopakumar-Vafa 
expansion plus non-perturbative terms, with respect to the string coupling parameter. 
It is here where we establish a remarkable exact coincidence with the conjectured non- 
perturbative completion of the topological string on the resolved conifold. In Section 
m we extend the results to orthogonal groups. Section |5] is devoted to the quantum 
limit of [19] applied to rehned Ghern-Simons, corresponding at large N to the well- 
known Nekrasov-Shatashvilli limit. In particular a S-dual like relation between the 
perturbative and non-perturbative parts of the free energies will be discussed. Some 
more general remarks will also be given in this section. In appendixes we present some 
series identities and supplemental details about multiple gamma and sine functions, 
used heavily in the body of the paper. 
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2 Universal Chern-Simons 


2.1 Generalities 


Recall that the universal Chern-Simons free energy reads [TO 


with 


and 


7 ^ := 


dx 


x(e* — 1) 


7 = r - 7 


- d) 


■II 


^\og{t/6) , 


( 2 . 1 ) 


7^^ : = 


dx 


x{e^ — 1) 


{f{x/t) - d) , 


f{x) := 


d := 


sinh 


x{vi —2t) 
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sinh 


x{v2—2t) 


sinh 


x{v3-2t) 


sinh (2f) sinh (2f) sinh (2f) 
(ui - 2t){v2 - 2t){vs - 2t) 


( 2 . 2 ) 


ViV2V^ 

where we dehned the effective coupling constant 5 := k + t with k the usual Chern- 
Simons coupling constant. The parameters Vi occurring in f{x) and d are referred to 
as Vogel’s parameter. For particular choices of Vi one can recover from fl2.ip the free 
energy of Chern-Simons theory on with all simple Lie gauge groups. In particular, 
this requires to impose Vogel’s condition t = Vi + V 2 + v^. Then t is identified, in a 
special normalization, with the dual Coxeter number h of the corresponding simple 
Lie algebra. Normalization mentioned is called minimal one and is defined by the only 
negative Vogel parameter (usually Ui below) to be equal to —2. As discovered in [T3] . 
m also includes the rehned Chern-Simons theories of mi at appropriate values of 
parameters, though Vogel’s condition will not be satished anymore in the rehned case. 

It is convenient to rewrite 7 as follows. We redehne x ^ txj 6 m 7^^ such that 

^ dx 


7^^ = 


Using the relation 


— 1 ) 


^ax _ gbx 


Ui.x/6) - d) . 


sinh 


x{a—h) 


ohx 


(2.3) 


g.. (e“^ - l)(e''"^ - 1) 2sinh (f) sinh (f) ’ 

and making use of that the combined integrand is even under x —x, we can write 

7 as. 


^ =-flog (i/'s)+ i 


dx sinh (a;(t — h)) 

X sinh {xt) sinh {x6) 


{f{2x) - d) , 


(2.4) 
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where we deformed the integration range to pass the origin on an inhnitesimal semi¬ 
circle in the upper half of the complex plane. We will refer to the deformed contour as 
M+. (We could have equally well deformed to pass along in the lower half, which we 
will denote as M_). We further rescaled x —)■ 2x5. 

As we will show in section 12.21 the integration over the dy term in fact cancels 
against the log term in fl2.4p . such that we obtain the (novel) neat expression 


1 r dx smh{x(t-S)) j(Ox) 

4 X sinh {xt) sinh (xff) ’ 


(2.5) 


for the universal Chern-Simons free energy. One should note that the proof of equiva¬ 
lence between the integral representation (12.Ih . and so (12.bh . and the ordinary (refined) 
Chern-Simons partition functions, performed in m and na, is completely analytical 
and exact. Besides that it simultaneously encodes the usual and refined Chern-Simons 
theories in an unihed way, other benehts of the integral representation fl2.5p are that it 
naturally extends Chern-Simons theory to non-integer and/or negative values of N and 
the refinement parameter f3, to non-integer values of 6 (c/., [20]), and more generally 
to wide ranges of complex values of parameters. Particularly it allows an exact large 
N expansion without any need to perform a semi-classical analysis. This, in particular 
the last point, will become more clear below, where we will discuss ways to explicitly 
evaluate the integral occurring in (12.5p as a sum of residue. 

Note hrst that it is not easy to establish convergence restrictions for the integral in 
fl2.5p for arbitrary values of parameters. Here, we make the assumption that 2t > Wj, Vi 
and t > 0. In particular, this assumption is satished for the parameters corresponding 
to SU{N) and SO{N) (with (3 > 0). We further restrict to k > 0 corresponding 
to 6 > t > 0 (we can also discuss k < 0 similarly, but omit it here). Denoting the 
integrand of fl2.5p as X, we can estimate under the above restrictions on parameters 
that for large x one has X ~ iglThis leads us to the condition (for the 
real parts) 
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K > 0 : 2t — ^^(d + |d|) < 0 

i=l 


( 2 . 6 ) 


needed for a convergent integral. For parameters of rehned SU{N) and SO{N) this 
condition is indeed satished, c/., [13] and the example sections later on. 

In order to show that one can directly rewrite the integral fl2.5p as a sum of residue. 
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(2.7) 


T = ‘I'Ki Res(X, , 


where {x^} denotes the set of poles of the integrand, one has to establish that one can 
deform the integration contonr withont picknp np an extra contribntion. That this 
is indeed the case follows from the mnltiple sine representation we will discnss below. 
One shonld however note that the pole structnre of the integrand is snbtle. Depending 
on the particnlar valnes the parameters take (e.^., integer mnltiplies of each other), 
enhancement to higher order poles may occnr for some snbsets of poles. 

2.2 Multiple sine representation 

Thongh fl2.7p looks at hrst glance very simple, the expressions resnlting from taking 
residne directly of X are in fact not as easy as one might expect. In particular, one 
needs to invoke non-trivial summation identities to bring the result to canonical forms, 
used in literature. Following ideas of [12] , it is however possible to hnd a closed simple 
expression for (12.hh in terms of multiple sine functions for which the integral represen¬ 
tation, and so the residue calculations, are somewhat different (but of course equivalent 
due to non-trivial identities), but already more or less in canonical form. The price to 
pay is the introduction of some mathematical machinery of multiple Barnes’ gamma 
functions, multiple sine functions and their integral representations. This, however is 
justihed by compactness of the resulting expressions, clear rules of transformations and 
mathematical rigorousness, in particular justifying that simple residue taking of X is 
indeed valid. 

Definitions 

We need hrst to recall some basic dehnitions. Barnes’ multiple (^-function IZH is dehned 
as 


r>r\ 



Note that Cr(' 2 ,s|;u;) is well dehned if all parameters Wi he on the same side of some 
straight line through the origin, and Re^; > 0, Res > r. In similarity with the dehnition 
of Euler’s gamma-function in terms of Riemann’s zeta-function, one can use Cr(^; s\w) 
to dehne a multiple gamma-function 


lnW(z|;u;) := •••, tf^r) |s=o • 
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This definition follows [22] and differs from the original Barnes’ one [21] by some 
modular “constant”, depending on parameters. 

Next, the multiple sine functions are defined via Barnes’ gamma function as 



— 2; 


r^(z| 

hi) 


(2.9) 

where |w| := Some useful identities for the multiple sine functions we will 

make us of in this section are 

Sr{cz\cw) = Sr{z\w) , ( 2 - 10 ) 

S'2(l|l,tc) = y/w. 


(More identities are listed in flB.7|) .i 

According to [23] [2l] the multiple sine functions Sr possess for r > 2, 0 < ReWj, 
and 0 < Rez < |;u;| an integral representation over the entire real line, bypassing the 
singularity at zero either in the upper or lower half of the complex plane, he.. 


m 


*S'r( 2 :|w) = exp ( { — 1Y — Br^r{z\’^ + (“1) 


dx 


T\ 


= exp 


TTl 


{-ly-^-ErM^) + (- 1 )^ 


dx 


1 ) 


^ nLi( 


oUJiX 


1) 


( 2 . 11 ) 


where 

function 


refer to the generalized Bernoulli polynomials, defined via the generating 


n 




1 ) 


'^^Br,n{z\w) . 

n\ 


72=0 


( 2 . 12 ) 


Note that the Bernoulli terms occurring in fl2.1ip are izvr times the residue of the inte¬ 
grand at x = 0, he., they are equal to half of the integral with same integrand over the 
small circle around x = 0. These terms are necessary to have an equality between the 
two different integral representation above. The mentioned integral representation of 
these Bernoulli terms leads to the similar to fl2.1Up identities for Bernoulli polynomials. 
In particular, we will need in the following discussion the identity [23] 


Br,n{z\w) = {-lyBr,n{\w\ - z\w) . (2.13) 

Using results of [23], the logarithm of multiple sine functions can be expressed as a 
sum of residue, he., 

logS'r(2;|w) = (—l)”^Rr,r(^|fh) + (—1)'’ Bes{S{z\w),xy . (2.14) 

Uhl 
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According to [23], this expression is valid provided Imz > Ini|;u;| > 0. (Similar for the 
contour R_, for which one needs 0 < Im^r < Im|w|). However, there is one subtlety. 
Actually, the bound of [23] appears to be slightly too strong. Namely, we believe that 
one can relax the bound to >, via an improved estimation, performing for instance a 
similar discussion as in [25] (see their section 3.5). 


Quintuple sine free energy 

We now have everything at hand to rewrite fl2.5p in terms of multiple sine functions. 
Let us however first show that the log term in (12.41) indeed cancels against the integral 
over fl2.3p . We can rewrite (using the equalities in fl2.3lL 


dx sinh(x(f —5)) 1 

X sinh (xt) sinh {x6) 2 


dx 


^xt 


- l)(e'5^ - 1) 


X 


1 

2 


dx 


^Sx 


X - 1) ■ 


Note that we have due to the identity fl2.13p 


5) - Br,ri5\t, 5) = {I - (-1)^5^,^(5|t, (5) , 

which equals zero for r even. Hence, we can freely add this pair of Bernoulli polynomials 
to the expression above, thereby being able to rewrite 


1 r dx sinh (x(t — (5)) 1 / S 2 {t\t,S) 

4 Jr X sinh (xf) sinh (x(5) 2 ^ \^S'2((5|t, 5) 




where we made use of the multiple sine integral representation fl2.1ip and the identities 
given in (I2.10p . We conclude that the log term in (12.4p indeed cancels out. 

Finally, let us discuss how to express the universal Chern-Simons free energy fl2.5p 
entirely in terms of multiple sine functions. For that note that the integrand X is a 
product of five sinh factors. We can expand the sinh factors in terms of exponentials, 
yielding 


y"/. 


dx 


d^a + \\w\)x 


X (e2xi - 1) (62^-5 - 1) nLi - 1) ’ 


where {a} is the set of 16 tuples of length 4 of all possible sign combinations, he., 
a = {(To,..., (Ts) with (Tj = ±1, s((t) := HLo parity of the tuple (s(u) is similarly 

defined), 


= 


cro(f - <5) + X - 2t ), 


i=l 


(2.15) 
























and Iwl = 2{t + 5) + Y^^=i \'^i\ ■ Now note that the set a has a Z 2 symmetry identifying 
tuples differing by overall sign. Hence we can write 


J^ = s{v) 

{<t}/Z2 


2 


g( 2 I —^ -j- ^(2! —I ^ 

X (^(,2xt _ X) (^e2x5 __ X) (^ex\vi\ _ X) ■ 


Clearly, the identities fl2.1ip and fl2.13p can be shifted hy z ^ z + ^\w\. In particular, 
we have 


Br,r{\w\/2 - Z^\w) + Brx{\w\/2 + = (1 + {-lY)Brx{\w\/2 - Z^w) , 

vanishing for r odd. Hence, we can again freely add such pairs of Bernoulli polynomials 
to the above summation. Making use of fl2.1ip and fl2.13p . we conclude 

J' = s(n)log JJ 5*5 + |w|/2| 2f,25, |ni|, • (2.16) 

{a}/Z2 

Note that the condition 0 < Re|z| < |xn| needed to rewrite the integral representation 
in terms of multiple sine functions is equivalent to fl2.6p . as is evident from fl2.15p . 
Denoting as above the integrand in fl2.1ip as iS(z|xn), we hnally arrive at 


17) 

where as before denotes the set of poles on the upper imaginary axis. Equation 
fl2.17p constitutes the main result of this section. The expression fl2.17p is simpler than 
(12.71) . because the numerator of the integrand 5(2;|xn) is not a product of trigonometric 
functions as in (12.7p . but just an exponential. In particular, due to the used multiple 
sine identities, the exponentials do not combine anymore to trigonometric functions. 
This implies the existence of non-trivial summation identities to map to the results 
obtained via the direct integration (c/., appendix lAP . 

3 Example 1: Ajy-i 

The hrst example we consider is rehned Ajv-i, whose universal Chern-Simons repre¬ 
sentation has been derived in [13]. Namely, the rehned Ajv-i theory sits at 

vi = -2, V2 = 2Y , V3 = t = YN . 


B = -s{v) ^ s{a) {Res {S{\w\/2 + zYw),x*^) +'^B^ Aim]/"^ + Za\vi, 

{xX},{a}/Z 2 ^ 


( 2 . 
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In order to compare with topological string theory, it is convenient to introduce the 
parameter 



PN 

~T 


(3.1) 


Usually, this is the t’Hooft parameter kept hxed at large N. Since we do not perform 
a large N expansion, for us n is just a parameter. The integrand of (I2.7jl specializes 
under this choice of parameters to 


^ sinh — 1)) sinh + 1)) sinh (a:(5/i — (3)) 

^ 4x sinh (S/ux) sinh (Sx) sinh (/3x) sinh (x) 

As for these choice of Vogel’s parameters one pair of sinh cancelled out, it is more 
convenient to write the corresponding partition function Z := in terms of quadruple 
sine functions, rather than via the general expression fl2.16p consisting of quintuple 
sines. Making similar considerations as in section 12.21 we almost immediately deduce 
(after writing the sinh in terms of exponentials) 


= Si{P\w) ^4(<5|w) 84(1 + (3 + (5|w) ^ 4(1 + 2(5/i|w), (3.2) 


with w = {1, (3,6,6fi). Using recurrent relations and identities of multiple sine functions 
(c/., appendix iBll. the partition function can be simplihed to 




1 Ssi^l + (jp|l, f3, 5) 

SMhf3.6) 


(3.3) 


Note that the multiple sine representation is convenient for investigating symmetries 
of the partition function. For example, at = 1 we clearly see the level-rank duality 
k N (up to the prefactor). Another potentially interesting symmetry is the exchange 
(3 kk 6 , provided t = 6 fi = (3N is invariant, he., if simultaneously N —)■ N(3/ 6 . 

The for us important symmetry of the partition function is with respect to the 
transformation 


(3^1/13, N ^ (3N + 1- (3, 6^6/13, p ^ p + . (3.4) 

0 

Since this symmetry involves inversion of (3, it turns via the relation between (3 and 
the equivariant parameters of the hl-background [26] 


11 = \/Pg,, (2 


1 


(3.5) 


into the hl-background symmetry ei -H- — 62 - 
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Generally, at large N the parameter p is viewed as fixed t’Hooft coupling constant, 
i.e., it is assumed to not transform under fl3.4l) . However, the /^-inversion symmetry is 
still manifest. Namely, shifting and rescaling 


h ^ h “ ^ (1 “ /^) ) 
5 —)■ 5 , 


(3.6) 


we have in terms of fixed /i and 6 

f3 


^A(h;/5) = 


Ssiliy/P + v(d, S) 


p- Ul- (3) 




(3.7) 


Up to a logarithmic term originating from the first factor, we conclude that the exact 
free energy 3Fa{.J 3', f3) is invariant under 

/?-)■-. (3.8) 

Note that the transformation (13.61) is the usual shift one has to perform to obtain an 
even powers of gs only expansion of refined free energies, cf., [271 EH]. 

Let us continue to evaluate fl3.3p . For that, we make use of the integral representa¬ 
tion fl2.1ip . yielding 


3Fa ~ 


X - 1) - 1) (e^^ - 1) ’ 


(3.9) 


where we dropped for convenience the log /i term and two generalized Bernoulli polyno¬ 
mials. We can solve the integral via summing over residue of the integrand, c/., fl2.14p . 
However, care has to be taken about what values the parameters take, as enhancement 
from simple poles to higher order poles may occur. In particular, for simplicity we will 
assume that 6 ^ Q. Note also, that as one ray of poles depend on 1/5 ~ there is a 
natural split into perturbative and non-perturbative poles, the latter being independent 
of 5. Hence, there is a natural split as well of the free energy, i.e., 




The residue calculation of fl3.9p is straight-forward, except that one has to distinguish 
in the non-perturbative sector between f3 rational or not. The former will introduce 
some additional calculational complications due to pole enhancement of some of the 
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poles. (Nevertheless the free energy for rational f3 can be derived as well by a non-trivial 
limit from the non-rational jS free energy.) 

Instead of stating directly here the resulting expressions for and it is 

instructive to calculate these free energies as well via the direct integration of Xa 
described in section 12.11 Though this approach is on a technical level more complicated 
than using multiple sine functions, it is to some extent more illustrative, as for instance 
flop invariance is explicitly manifest, the Bernoulli polynomials are absent, and we can 
easily invoke trigonometric identities. 

As we assume 5 ^ Q, so we have p ^ Q. Equation fl2.7p tells us that 

Ta = 27rz Res(XA, x*^_). 

{<} 

The set of poles consists of two rays of perturbative simple poles at x* i = rnii/^d^) and 
Xp 2 = rnii/d (since p ^ Z) and two rays of non-perturbative poles at x*^^ ^ = mii/{3 
and 2 = nni, where n G N \ {0}. Hence, 

P) + -^P,2 ; P) •= Rnp,l + Rnp,2 , (3.10) 

where R denotes the residue sum of the respective ray of poles. (The sums here and 
below are actually understood as a limit of sum of partial sums of all objects involved, 
since infinities are cancelled in between them.) The perturbative contribution Jyf can 
be easily inferred, up to some technical subtlety. Namely, only the combined RpA + Rp ,2 
is finite over summation over the set of poles. In particular, taking partial sums implies 
that we pick up a left over contribution of — | log fi in canceling the singularities against 
each other, which in fact exactly reproduces the prefactor of fl3.3p . Using trigonometric 
product-to-sum formula and addition identities, it is not hard to infer that we can write 
the combined residue summation as 


1 1 

V ^ ^ n=l 

cos (n7r(2/i + (1 — fi)/S)) 


nsin(^) sin(^) 

4 n 

n=l 

(3.11 


Note that after shifting the free energy via fl3.6l) . the /^-inversion symmetry is clearly 
visible, and furthermore the free energy becomes explicitly flop invariant {i.e., invariant 
under fl —fi)- 


In order to make contact with the free energy of the rehned topological string, we 














and make use of the (non-trivial) sum identities flA.3j) and flA.Sj) . such that we arrive 
at 


_ 0 ^ _ 1 g-ng.(V?-l/V?)/2 

i ^ nsinh sinh ^ 2 ^ ’ 

(3.12) 

where we also dehned Q := and dropped some terms consisting of generalized 

Bernoulli polynomials and as well the log jj, term. We recognize the rehned Gopakumar- 
Vafa expansion and constant map contribution of the resolved conifold [5] . 

Let us move on to the non-perturbative contribution. Note that we have to dis¬ 
tinguish between /3 G Q* and (3 ^ Q*, with Q* := Q \ {0}, as already mentioned 
above. 


If 13 is not a rational number, we just have two rays of non-perturbative simple poles 
at i = mri/(3 and x*p 2 = leading to 



n=l 


sin sin j 

nsin sin 

sin — 1)) sin — (3/6)) 

n sin {n7rf3) sin (nirS) 


(3.13) 


Via invoking trigonometric identities, and under the redehnitions fl3.6l) . we arrive at 


j-r(h;/3^Q)=TE 


°° ^ cos I ^ ^ + nTTj cos (?7.7rv^5(2/i — 1) + utt) 


=1 \ n sm I sin sin {n7r/3) sin {mr^S) 

, l + cot(^)cotf^ 


-lE 


^ J V {n7r(3) cot {mry/pS^ — 1 


n=l 


n 


n 


(3.14) 

where we made use of the periodicity cos((r + nvr) = cos((r — nvr). The inversion 
symmetry fl3.8l) is clearly manifest in the non-perturbative part of the free energy for 
/9 ^ Z, as expected. 
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It remains to write fl3.14p in terms of gg- Making use of flA.Sp and (lA.Sp . we deduce 






— 2n7r^ 27rm 

e Q Vpga 


AE 


— 2n7r‘^y/^ 27rm-\/^ 

e 9s Q 9^ 


AE 


n=i n sin j sink ^ n=i n sin {mij3) sink j 

off; 9 a) 


2“ nsin(^- 


(lA) 


AE 




/ 2n-w'^y/p\ 


n=i n sin irniB) sink . , 

^ ' \ 9 s J 


(3.15) 


(We dropped again some generalized Bernoulli polynomials.) 

/3eZ* 

As for our purposes, i.e., comparison witk known non-perturbative results at /3 = 1, 
it is sufficient to take /3 G Z*, we restrict for simplicity kere to tkis case. We kave 
one ray of non-perturbative double poles at x*p 2 = ray of simple poles at 

^np,i = rmri/(3 with m G M \ {1(3}. We infer 


= -) E 


sm 


n7n5(/i— 1 ) 
0 


sm 


n7ri5(^+l/5) 


sin 


CX) 

Rnp,2 ^ ^ 

n=l 

CO 


E 

n=l 

00 

E 

n=l 


n sm 

{1 — (3 + 25) — {1 — (3) cos (2n7r5) 

8n/3 sin^ (mrS) 

((3 — 1 — 25/2) cos {2mr5{/2 — 1)) — (/? — 1 + 25(1 — /i)) cos {2mr5/2) 

8n(3 sin^ {mr5) 

sin {2mi5) + sin {2mi5{/2 — 1)) — sin {2mi5/2) 


8n'^7i(3 sin^ {mi5) 


(3.16) 


Note that for (3 integer some of the simple poles are enhanced to double poles and 
therefore the above residue results are more complicated than before. In order to 
simplify Rnp ,2 further, we redefine /i /i'+ 1/2 such that with the help of trigonometric 
sum-to-product identities we are led to 


Rnp,2i/9‘ ) — ^ ^ 


n=l 


1 — (3 sin {2mi5) 

4:(3n Sn'^TiP sin^ {mi5) ' AnP sin^ {mi5) 


+ 


E 

n=l 


cos {2mr5/2') (P — 1 — 25/2') sin {2mr5/2') 5 cos {2mr5/2') cot {mr5) 


AtPttP sin {n7r5) 


+ 


AnP sin {mT5) 


AnP sin (utt^) 


(3.17) 
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Adding up the contributions and with some further trigonometric simplihcations, we 
conclude that 


p p 4nsin nTTd 

n=l ^ ' 


(3 sin (2mi6jjL') 


E 


cos 


n7n5(2^'+l/5) A 


P 




1 ^ cot cot 


n^i/3 n sm 


f) (t) 


-zE 


n^l(3 


n 


(3.18) 


with 

CX) 

z^tVa)=E 


n=l 


cot (n7r(5) 

' 4nsin^ (riTT^) 


+ 


/ cos {2nn6fi') ^ 6fi' 

f ^ \ cin I'n'TrAi 9' 


^ sin (riTrh) 2?7,sin(n7r5) 


sin {2miSfi') ^ 6 cos {2mi5^') cot {mi5) 


An sin {mi5) 


(3.19) 


The /5-dependent non-perturbative completion for /3 G Z* given in 03.181) appears to 
be non-trivial and constitutes one of the main results of this work. Note that the 
symmetry 03.8p is not visible, as it maps /5 gZ*—)-1//5gQ. The symmetry however 
dictates the solution for 1/(3 G Q) with (3 ^ Z*. 

Finally, let us inspect 1) in some more detail. The exact resolved conifold 

free energy at /5 = 1 has been calculated already previously in [16]. However, it is 
interesting to compare the exact result predicted by universal Chern-Simons to other 
recent predictions. For that, note that under usage of the identities in section IATI and 
in terms of Qg, we have (up to some generalized Bernoulli polynomials) 


1) ~ E 


n= 

CO 

E 


t 47^^r^2 sinh 


t 47rzn2 sinh 


1 + 


1 + 


ry 2n7T^ i 

2 ?7,7r^ e 3s 

AnTr’^ifi' 2mr'^i 


e 3s 


coth 


2n'K‘ 


ijj, 

e ^ 


(3.20) 


This has to be compared with the prediction for the non-perturbative completion of 
the topological string on the resolved conifold inferable from the conjecture of HI- 
Under setting 
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we precisely recognize in the second row of fl3.20l) the conjectured M2-brane contribu¬ 
tion, stated for instance in [18] (c/., their eq. (5.40)). The first row corresponds to the 
non-perturbative completion of the constant map contribution, which, to our knowl¬ 
edge, is not given explicitly in the literature. We conclude that the exact calculations 
via universal Chern-Simons give strong evidence about the validity of the conjecture 
of [T7| for the resolved conifold. Turned around, this implies that the large N duality 
is exact. (Very recently, while our work was in the writeup process, some numerical 
checks of the exactness of the large N duality between Chern-Simons and the unrefined 
conifold appeared, see [29].) 


4 Example 2: D]\i /2 

In this section we will discuss the free energy of refined Dn /2 Chern-Simons theory. 
Invoking the multiple sines technique introduced in section 12.21 the free energy can 
be obtained via a simple residue calculation as a Gopakumar-Vafa sum, and non- 
perturbative corrections thereto, similar as in the previous section. 

The parameters leading to SO{N) refined Chern-Simons theory {N is understood 
to be even) have been derived in [13] and read 

Vi = -2, ^2 = 4 / 3 , V 3 = P{N-4), t = P{N-2), 6 = K + t. ( 4 . 1 ) 


The multiple sine expression of the free energy T derived in section |2^ i.e., fl2.16p . can 
be used for calculation of free energy in this case. However, without passing through 
quintuple sine functions, we can infer via a more direct independent calculation that 


^3(l + |t|l,t,5) S,{l + ^ + t)\l,2P,6) 
^3(l + 5+|t|l,C5) "" SMl,2P,6) 


(4.2) 


where Zd{N] (3) := e . The first factor can be simplihed as 


83(1 + _ 82(1 + _ S'2(|t|l,t) _ ^/2 

^3(l + <5 + it|l,C5) “ ^2(l + |t|l,5) “ ^l(|t|t)^2(l + |t|l,<5) “ ^2(l + |t|l,<5) ’ 

where we made use of the identities 8i{t/2\t) = 1 and 82 {t/ 2 \l,t) = \/2. Note that 
the remaining non-trivial denominator S' 2 (l -|- |t|l, h) becomes S' 2 (|iV|l, 5) in the non- 
rehned limit (3 = 1 and matches one of the multipliers of the partition function of 
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SO{N) Chern-Simons given in [16] (recall the invariance of multiple sines under simul¬ 
taneous rescaling of the argument and all parameters, cf., fl2.10p '). he., 


7 IN u-o-i / g3(2iV|2,2,2.5) ^S,{2N\i,2S) 

' V S3(2|2,2,2i) S2(iV|2,2i) 


(4.3) 


Let us verify that the other terms in (14.31) appear from the remaining second factor 
in the rehned partition function fl4.2p . For that, note that the denominator of the 
second factor in fl4.2l) can be rewritten for = 1 as (see appendix [B]) 


^3(1|1,2,<5) = V^3(1|1,1,<5)V^2(1|2,<5) = 2V^3(l|l, 1,5), 


which exactly matches the corresponding term in the SO{N) Chern-Simons partition 
function given above. So, our rehned SO{N) partition function equals the SO{N) 
Chern-Simons partition function at /3 = 1, as it should be. 

Let us move on to the rehned case with general f3. The multiple sine representation 
(14.21) of the rehned partition function is essentially (we omit numerical multipliers) 


Zd{N;P) ~ 


^3(l + /g(iV-l)|l,2/3,^) 

^3(/5|1, 2(3,6) ^ 2(1 - i/9 + - 1)|1, <5) ■ 


(4.4) 


This form of the partition function is most suited to invoke the integral represen¬ 
tation of multiple sine functions (12.lip in order to obtain the free energy as a sum 
of residue {cf., fl2.17p ). However, we do not need the full expression fl4.4p for that. 
Namely, if we dehne 

I3N -(3 

h:- , 

and compare with Za given in (13.3p . we infer that actually 


Znip-, P) ~ ZAm 2^^) X M(P) 


(4.5) 


with 


T{^^■,P) := 
M(/3) : = 


^log/i - logS'2(l - (3/2 + (5/i/2|l,5), 

^3(2/911,2/3, ^) 

^3(/9|1,2/3,<5) ’ 


(4.6) 


holds. (We made use of the identities (lB.7p ). Note that /i is the usual t’Hooft coupling 
kept hxed in a large N limit. However, we like to stress that we actually do not 
take any limit. We observe that the rehned Dj^/^ partition function splits into an 
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oriented and unoriented sector, c/., [30]. Usually, T is referred to as domainwall tension. 
The function M{f3) is the unoriented piece of the constant map contribution and is 
expected to relate at /5 = 1 to the real MacMohan function introduced in the context 
of orientifold in [3ll [32] . 

As a side remark, note that the hl-background symmetry discussed in the Ajy_i 
example does not carry over to Dj^f/ 2 , as the domainwall tension (or orientifold plane) 
breaks the symmetry. Instead, it seems that the related transformation 

1 / 4 / 3 , 6 ^ 6 / 2 / 3 , 


combined with some appropriate transformation of N, maps partition functions of 
rehned SO{N) to rehned Sp{N) and vice versa. Since we do not have at hand the 
rehned universal version of the Sp{N) Chern-Simons partition function, we can not 
explicitly verify this hypothesis at the time being. (However, one may use the expected 
transformation properties, and the expected structure of the rehned partition functions, 
c/., [Sn], to predict the corresponding Vogel parameters.) 

As we already discussed the non-perturbative completion of the Gopakumar-Vafa 
expansion of Z^(yU,;/3) in the previous section, we only need to discuss here the addi¬ 
tional contribution due to the factors M and T. Let us start with T. From the integral 
representation (I2.1ip we immediately deduce 


1 ? TT 

^ = 2 + ^-62,2(1 — / 3 / 2 . + (5/i/2|l, 6 ) + 


gx(l-/3/2+<5M/2) 

X (e^ — l)(e'^* — 1) 


(For convenience, we drop from now on the hrst two terms.) As in the previous Aj^ 
discussion, we can distinguish between perturbative and non-perturbative poles (the 
former depending on Qs). Hence, there is a natural split into a perturbative and non- 
perturbative part, he., 

r = r^ + T^p. 


As we assume 6 not to be integer, there are however no cases to be distinguished. 
Hence, we can write down immediately 


r(/x;/3) 




E 


GE 


2 -^ nsin(^) 2“ nsm(mr6) 

n=l \ 0 / n=l ^ ' 


(4.7) 


where the hrst sum corresponds to and the second to T 


■NP 
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For later reference, note that T can be as well expressed as a linear combination of 
the free energies Ta- This can be seen as follows. Using the second identity listed in 


flB.7ll we infer 

r(M;/3) ^ ^3(l-/g/2 + & + W2|l,&,<5) ^ Z^(/i/2 + (46-/3)/(2^);2&) 

"" ^3(l-/5/2 + 5/i/2|l,^<5) ZA{fi/2-p/{26)-2b) ’ 

where 6 7 ^ 0 is arbitrary. Hence, taking for instance b = 1/2, we conclude 

r(/i; f3) = J^Aib^n + (2 - / 3 )/( 2 < 5 ); 1 ) - J-^(/i /2 - / 3 /( 2 < 5 ); 1 ). (4.8) 


It remains to discuss M{I3). Ignoring for notational convenience the generalized 
Bernoulli polynomials, we have 

f diX P dx 6 ^^ 

logM(/3) ~ + • 

As is by now familiar, we have a split into a perturbative and non-perturbative piece. 
Summing the residue of the perturbative poles, yields 

. 00 inn{P-l) 


logM'’(/J) = jJ] 


4^nsm(f)cos(:f)' 

For the non-perturbative parts we have to distinguish between /3 G Q or not. 


We obtain 





gn7ri(/3-l) 

nsin (riTT/S) cos {rmS) 


1 

2 


E — 

nodd nsm 


i7m(5+l) 

e 23 




Note the non-perturbative appearance of an odd sector in n. For /3 G Q let us just 
give one example. Namely /3 = 1 corresponding to (unrehned) topological string ori- 
entifolds. 


/3 = 1 
We infer 




imrS 

2 


n=l 


n sin (^) 
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Hence, combining all parts, we deduce via relation fl4.5p and under using fl3.18p that 


1) ~-j-r(h; 1) + 1) 



72=1 


sin {2mi5 /i') 
nsin(n7r5) 



72=1 


sin(n7r(5/i') 
n sin (^) ’ 


(4.9) 


where the last two terms are the remaining terms of fl3.18p at /5 = 2 (the fractional 
sums can be resummed for this value of (3). The hrst of them, can be eliminated via a 
further shift of /i'. 


5 Quantum limit 

It is interesting to consider the quantum limit introduced in [19] at hand of /3-ensembles. 
Namely, taking N large with 

Af:=PN, 

hxed. In order to keep Af hnite, this also requires taking /3 —)■ 0. Hence, we dehne 


oW := hm/3J^. (5.1) 

(The notation of quantum limit originates from the fact that in this limit /3-ensembles 
are captured by ordinary quantum mechanics. However, this notation may not be the 
best choice as it rather corresponds to a classical limit reducing a double quantized 
system to a single quantized system.) 


An 


Let us apply the quantum limit to the rehned An case discussed in section [3l However, 
it is instrutive to apply the limit not after residue taking, but directly to fl2.5p . We 
infer 


= —J- 


dx sinh — 5)) sinh {x{J\f -|- 1)) 
2 


/]R+ x^ sinh (a;5) sinh (x) 

We have one ray of perturbative simple poles at x* = mri/d and one ray of non- 
perturbative simple poles at x*^^ = mii. Taking residue leads to 
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with 


„WX(A/-,i) 



n=l 


cos j 

n? sin (^f) 

cos {mi{2N' — 5)) 
n? sin (mrS) 




cot (f) 




E 

n=l 

1 cot (mrS) 


n=l 




(5.2) 


(Strictly speaking, the discussion in section 12.21 about validity of residue taking does 
not directly apply to the /3 ^ 0 limit. However, it is not hard to convince ourselves 
that the discussion can be adjusted to deal with (15.11) . as in this limit the integrals are 
actually stronger suppressed.) 

It is interesting to observe that the perturbative and non-perturbative part of oW^ 
are related by 


■5) = J <sWa fw - 1, i) 


(5.3) 


We will come back to this S-dual like relation between the perturbative and non- 
perturbative part of the quantum free energy below. For now, let us take an additional 
large M limit, keeping 

/i = y , (5.4) 

hxed, similar as in [19]. We infer that 


oWA((j/x) — A,Ns{f^) ) 


(5.5) 


where we defined the Nekrasov-Shatashvilli (for short NS) free energy [33] (see also 
[34] ) as 

J^A,Ns{l^) := lim /3 J) 4 (/i; /3 ), 

13^0 

but there we take now, in contrast to before, the limit after residue taking. The limit 
can be easily applied to fl3.1ip and (13.141) . For that, note that the summation index of 
the summations in (13. Mh which can be redefined to run over n/(3 move off to infinity 
in this limit, and therefore these summations can be dropped. Hence, 




A,NS 


(f) 


-pNP 
A,NS 


(f) 


6 cos (n7r(2/i-|-1/(5)) 6 cot (^) 

2tt ^ n3sm(^) 2n ^ m? 

n=l \ d / n=l 


1 COS (n(5(2yU — 1)) 
27r 77,2 sin (nvr^) 

n=l ^ ' 



n=l 


cot (httS) 




(5.6) 
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Comparison with fl5.2p confirms fl5.5p . Note however that if we had applied the limit 
instead to (13.181) . which is only valid for (3 E Z, we would have gotten a different result. 
In particular, we learn that 1) ^ Instead, it is easy to verify that 

the relation 


1) = N 4 


1 KF: 


■NP 
A,NS 


( 4 ) 


(5.7) 


holds. This relation is an exact statement for the non-perturbative part of rehned An 
Chern-Simons on free energies at large N. If we believe in the exactness of the large 
N duality, this translates to a property of the non-perturbative (rehned) topological 
string free energy on the resolved conifold. Making use of (15.3p . we can as well express 
-^4^(4; 1) via J^a,ns leading to the conjecture of [T^ (see also [T8]h 


D 


N/2 


Let us ask what happens in the Dn /2 case. Applying the quantum limit dehned in 
(15.ip to the integral representation specialized to rehned Dn/ 2 , we observe that 


,Wd{M) = ]^oA^a{M) , (5.8) 

where we kept M = fiN — (3 hxed. In fact, this non-uniqueness property of the NS limit 
has been already observed in [35] and [30]. Hence, trivially the relation (15.3p holds as 
well for Dn/ 2 - However, this is only one half of the story. As the /^-inversion symmetry 
is broken in the Dn /2 case, we have to consider as well the alternative quantum limit 

oo>V ;= lim ij'. (5.9) 

/3^oo p 

We also dehne 

Af = N-2, 


in this limit. (Note that in contrast to the (3^0 case we keep N hnite). Physically, 
after taking J\f large with (15.4p hxed, the occurrence of two diherent limiting cases 
can be explained by the fact that the orientifold plane lives in two dimensions of the 
space-time. As the Nekrasov-Shatashvilli limit corresponds to a reduction to two space- 
time dimensions, we have two choices. Either we reduce to the subspace hlled by the 
orientifold plane, or to the orthogonal subspace. The latter leads to the relation (I5.8p . 
as we do not see the orientifold in space-time anymore. The other case, with orientifold 
plane, can be investigated as follows. 
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The limit fl5.9p can be calculated via redefining x ^ x/ (2/3), and rescaling 5 ^ 213 5 
in the integral representation fl2.5p (c/., [I3]), leading to 

If dx /sinh (a;A/'/2) sinh (x(A/'/2 — (5)) sinh. {x{J\f/2 — 5)) 

^ 2 J^+ x"^ \ sinh (x) sinh {5x) sinh {x6) 



We denote the integral over the last term as ooT = ooT^ + ooT^^, and evaluate via 
summation over residue 


.T^ = - 


r- CO 

S sm 


f mrAf^ 
\ 25 ) 


IT 


n=l 




oor^^ = o. 


i 6 



Li 2 



(5.10) 


Note that coT^^ vanishes because the second term has only one ray of (perturbative) 
poles. This is consistent with fl4.7p as in this equation either or can survive 
the quantum limit, depending on how we rescale 5 before taking the limit. Evaluating 
as well the integral over the first term, we conclude 




r °° „„„ (n-KM\ j. oo 

A V + — V ^ , n-P 

2t[ ^ m? sm. (Pf-') 2 t: ^ m? sm (Pf-') °° ’ 

n=l \ 5 / n=l \ 5 / 


1 COS {rniiN' — 
2t\ ^ -n? sin iimS) 

n=l ^ ' 




n=l 


cot (nTrd) 




(5.11) 


Note that both, the domainwall tension and the constant map part appear not to 
be compatible with a relation like fl5.3p . 

Finally, let us ask if we nevertheless can still find a relation like fl5.7p in the T) 7 v /2 
case. For that, recall from section H] the relation fl4.5p . Hence, we have for instance 


jT(/i; 1/2) = J-r (h; 1) + iogM^^(i/2) + . 


(J~NP does not vanish away from the Nekrasov-Shatashvilli limit, c/., fl4.7p .l Recall as 
well that T can be expressed as a linear combination of 1) with shifted n via 

fl4.8p . (This corresponds to the well-known fact that the orientifold contribution on the 
conifold is essentially given by half a closed period, see [36].) Since fl5.7p extends to 
shifts n h + f) we deduce that 1/2) can be expressed via a linear combination 

of first order linear differential operators in 5 acting on shifted Nekrasov-Shatashvilli 
free energies up to the constant map contribution. 
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Remarks 

Several remarks are in order. Though the relation (15.dh looks very appealing, it rather 
appears to be a special property of the quantum limit of the rehned topological string 
on the resolved conifold (and other genus zero geometries), as such a relation is not 
obviously visible for the non-perturbatively corrected /3-ensemble [19] and M = 2 
supersymmetric gauge theory 1371 (see also [3H1I39]) quantum free energies (where the 
non-perturbative instanton corrections are due to B-cycle instanton tunneling, at least 
at a particular point in moduli space). In particular, a general validity of the S-dual 
like relation (I5.3|l would trivialize the theory of resurgence. 

However, the quantization condition used in the context of toric Calabi-Yau man¬ 
ifolds in BDl, see in particular HU, seem to invoke besides the known perturbative 
part of the quantum periods [3l| only the S-dual piece of the non-perturbative part of 
the quantum (or NS) free energies. As this is directly linked with the validity of the 
conjecture of im, i.e., that the relation (15.7p (combined with (15.3p ) in general provides 
the non-perturbative completion for topological strings on (general) toric Calabi-Yau 
manifolds, for which various numerical checks seem to have been performed in the 
literature, a puzzle arises. Namely, how can this be ? 

Possible resolutions could be as follows. Either there are in fact only A-cycle type 
instanton corrections (c/., |l2]) on toric Calabi-Yau manifolds, see however [13]. (Some 
might see the fact that the poles in the Gopakumar-Vafa expansion arise only at tree- 
level, c/., [H], as evidence.) There is as well a hidden S-dual like relation between 
the perturbative and non-perturbative parts of the (non-trivial) quantum free energy 
calculated in [191 (and so in ordinary quantum mechanics, which would actually be 
too good to be true). Or, perhaps most appealing, that the dominant instanton action 
changes over moduli space, i.e., at large volume A-cycle type instantons are dominant, 
while for instance at conifold points B-cycle instantons are dominant, c/., [HI [371 1^ - 
In the sense that the above simple relations break down away from large volume, due 
to quantum effects. We believe clarifying these points will be of general interest. 
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A Trigonometric sum identities 


A.l cos and sin sums 


Consider the summation 


CX) 

S+(x|w) := 

n=l 


cos (27rnx) 

ni=i 


From the dehnition of the generalized Bernoulli polynomials fl2.12p we infer the relation 


-| / • W_ f 

= (A.1) 

Note that in order to arrive at flA.lD we also made use of the identity Br^i{cz\cw) = 
d‘~^Br^i{z\w) {cf., [23]). Hence, we can rewrite as 


CX) 


E 


{2^YBr,l{\w\/2\w){^y-^ 

n 


OO 

cos(27rna;). 

n=l 


Recall the polylogarithm Lis( 2 ;) := ^ (dehned via the convergent series for \z\ < 

1, but can be extended to the whole complex plane via analytic continuation), and 
Jonquiere’s inversion formula valid for n > 0 (and 0 < Rex < 1 if Imx > 0, else 
0 < Rex < 1) 

LR(e2--) + (-l)-Li„(e-2--) = , 

nl 

with Bn{x) the ordinary Bernoulli polynomials. The above inversion formula extends 
to n < 0 if we dehne that i?„<o(x) = 0. Hence, in terms of the polylogarithm we can 
write 




(2iY ^BrY\w\/2\w)(iy ^ 


1=0 


X 


((1 - (-1)"«-') LW._, 


(A.2) 


Now note that BrY\w\/2\vi) vanishes for I odd due to the identity fl2.13p . But if I is 
even only, we have 


S., = 


{x\w)+Er= 


OO {‘2irBr,i{\w\/2\w)ii)^-- 


(x\w) 


Lh+<j-z (e ^^“) r + s odd 
r + s even 
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where we defined 


S(x|;u;) 


OO 


E 

1=0 


^r+s-l22r+s-l-l^r+s-l^^^^ /2 \w)Br+s-l (x) 

l\{r + s-l)\ 


Note that the summation in H is finite (as we defined the Bernoulli polynomials to 
vanish for negative n). Writing the polylogarithm again as a series and using flA.ip . 
we conclude with the identities valid for r + s odd 


E 

71=1 


cos (27rnx) 


ri' 


n 


i=i 


sm [nwi) 


S(x|^) + Er=i 

^(1-^N +Er=i 


n® 01=1 sm{nwi) 

^ 2 TTinx 

n® 01=1 sm{nwi) 


(A.3) 


where the equality in the lower row is due to the fact that we could have equally 
expressed E_|_ in terms of Li„(e’''“) (making also use of the identity {— l )"' Bn { x ) = 
Bnil - x )). 

Similarly, we can consider 


CX ) 

E_(a:|w) := ^ 

n=l 


sin {27rnx) 

ni=i 


The only difference to E_|_ discussed above will be a flip of sign and an additional overall 
factor of — h Hence, we have for r + s even 


E 

71=1 


sin {27rnx) 

riLi 


— iE{x\w) 
z5(l — x\w) 


-*E 

+*E 


QQ g —27rma: 

71=1 ri'® sin(nit;/) 

QQ ^ 2 'Kinx 

n=l 71^ n[=i sm{ 7 iwi) 


(A.4) 


As a side remark, note that for the parity of r + s there the polylogarithms cancel out 
we provide via the function H a finite expression for some of the summation formula 
in [25] . 


A.2 cot sum 

Consider 


s-E 


1 + cot(?7,7rx) cot{miy) 


n 


E 


1 — cot (nvrx) cot {—miy) 


n 


n=l n=l 

For simplicity, we assume in the following that Re?/ = 0. Invoking the identity 

2T[i 


TT cot TIZ = in -\- 


o 2 'Kiz _ 2 ’ 


(A.5) 
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leads to 




2e 


—2n7zix 


^^2n7Tix _ _ ^—2n'Kiy'^ _ ^—2n'Kix'^ _ ^—2n'Kiy'^ 

For Im?/ > 0 we infer with help of the geometric series 


S = 2 E 

n=l 

oo 

E 


32n7r2X 


- 2 E 


E— 

< ^ 71 C:iT 


^-mri(x-y) 


= I 


^ I g^,>„ ^ ^^ 2 nmy _ ^ gjj^ (jmx) sin (mry) 

sm{mr{x — y)) ^-n-Ki{x-y) 


V — 

( ^ 71 I51T 


^ n sin {mix) sin {mry) n sin (riTra:) sin {mry) ’ 

where we made use of fl2.3p . Similarly, for Imy < 0 we have 


E = z V — 

f ^ 71 151 


sin {mi{x + ||/|)) 


V — 

< ^ 71 I51T 


^-n-Ki{x+\y\) 


“ nsin {mix) sin {mi\y\) nsin (nvrx) sin {mi\y\) 

Invoking the summation identity flA.4D . we finally deduce (for Key = 0) 

+ Im|/>0 

2 nsAZ^hy, + S(x + \y\\x, \y\) lmy<0 


(A.6) 


(A.7) 


(A.8) 


B More on multiple gamma and sine functions 

Barnes’ multiple zeta functions (as defined in fl2.8l) ). fulfill recurrent relations pTl |22] . 
which were implicitly widely used in this paper. Indeed, for z = zq + t(;i,Re2;o > 0 
the sum over n* effectively starts from Uj = 1, with zeta function argument being zo- 
Adding and subtracting the contribution of n* = 0, we get the relation 


CN{Zo + Wi,s\w) = Cv(^ 0 , s|w) -CN-l{Zo,s\Wi,...,Wi_i,Wi+i...,WN), 


(B.l) 


In turn, from this follows the recurrence relation on multiple gamma-functions: 

Tn{zo\wi,W2, ...,Wn) 


Tn{zo + Wi\Wi,W2, ...,Wn) = 


r7V-l(2;o|'«^l, •••, UJi-l,UJi+l---,UJN) ’ 


Another type of relations between multiple gamma functions appear when there 
are some relations between parameters. This is based on the integral representation 
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[22] (see fl2.12p for the definition of generalized Bernoulli polynomials) 


log T]^{Z\w) 




1 

(1 - 


N-l 



n=0 



BN,n{z) 


(-1)^ 

N\ 





(B.2) 


and the observation of m that if in the linear combination of logarithms of some 
multiple gamma functions the main terms (he., the hrst terms in flB.2p ) cancel, then 
all other terms cancel as well. This yields a relation between the corresponding multiple 
gamma functions. As an example, consider the identity 


1 / 1 1 \ __ 1 _ 

2 \(1 — e~^y 1 — (1 — e“^)(l — 6“^*) ’ 

which implies identities between multiple gamma functions like 


r2(iv|i,2) = v/r2(iv|i,i)ri(iv|2), 

r3(iV|l,2,a;) = v/r3(iV|l,l,x)r2(iV|2,x), 


etc.. 

Evidently, such relations can appear each time one has rational relations between 
parameters. We also note that from the identity [22] 


ri(tc|a) = exp 





it follows _ 

ri( 2 :|x) = ri{x\2x) = ^, ro{w) = l/w. 

where the last identity holds by dehnition. 

Recalling the dehnition of the multiple sine functions in terms of gamma functions 
fl2.9p . it is now easy to show that we have 


^3(iV|l,2,<5) 


1 _ r2((5 + 2-Ar|2,5) 

r 3 (iV|l, 2, <5)r3(<5 + 3 - iv|i, 2,5) “ r 3 (iV|l, 2 , 5)Ts{5 + 2 - iV|l, 2,5) 

_ T2{S + 2-N\2,6) _ 

Vr 3 (iv|i, 2 , <5)r2(iv|2,5) v/r3(<5 + 2 - iv|i, 2 , <5)r2(<5 + 2 - iv| 2 , 5 ) 
V^3(iV|l,l,<5)V^2(iV|2,<5). 


(B.4) 





















This identity is made use of in section HI 

Clearly, as similar integral representations as flB.2|) exists for the multiple sine func¬ 
tions (c/., (12.lip ) [23], one can as well derive new identities for multiple sine functions 
in a similar spirit. For instance, from the relation 


+ 


(e^ - 1)2 e 

we can deduce the identity 


2x 


- 1 


- l)(e2^ - 1) ’ 


(B.5) 


S2{N%y) 


= S^{N + l\l,2,y). 


(B.6) 


As a side remark, note that for the multiple sines it is even more evident that relations 
of the type flB.Sp lead to identities between multiple sine functions, since both terms 
in the integral representation, he., integral over real line and generalized Bernoulli 
polynomial (the residue contribution of the origin), are integrals of the same integrand 
of type fIB.Sp . 

Finally, for the readers convenience, we list some of the simpler identities of sine 
functions we made heavily use of in the main text 


Sr{cz\cw) = Sr{z\u) , 


Sr{z + Ui\u) = Sr{z\(^ / Sr-l{z\!^^ ), 
Sr{z\u) = S'r(|a;| — , 

S'i(z|a;) = 2 sin — , 

UJ 

S 2 {l\l,x) = y/x, 


were := (cui,..., cUj.i, cUi+i, ...cu^). 


(B.7) 
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